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Abstract 

We report exact multielectronic ground-states dependent on electron concen- 
tration for quantum mechanical two-dimensional disordered two-band type 
many body models in the presence of disordered hoppings and disordered 
repulsive finite Hubbard interactions, in fixed lattice topology considered 
provided by Bravais lattices. The obtained ground-states loose their eigen- 
function character for independent electron approximation, perturbatively 
are not connected to the non-interacting but disordered case, and describe 
a localization-delocalization transition driven by the electron concentration, 
being highly degenerated and paramagnetic. 
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I. INTRODUCTION 



In the real life, the crystallin state is the exception rather than the rule 1 , and as a 
consequence, the disorder exists everywhere, ranging from few impurities or interstitials in 
a periodic host, up to the completely disordered glassy and amorphous structures, alloys 
and compounds. Given by this, the effects of the disorder are intensively analyzed 2 , special 
attention being given in the last period to two-dimensional (2D) systems, where the ob- 
servation of metallic behavior in 2D high-mobility samples 3 contradicting the conventional 
non-interacting scaling theory 4 , has underlined the special importance of electron-electron 
interactions in disordered systems, at least when its value is relatively high 5 , or when the 
competition between disorder and interaction demands the consideration of both 6-12 . 

Deep-rooted in the difficulty of describing the effects of the disorder in a non- 
approximated manner, on the theoretical side the interpretations are given almost exclu- 
sively based on approximations. In the last decade however becomes to be clear that this is 
not a fortunate situation, since not only the non-interacting scaling theory has been affected 
by new experimental results, but also other approximated schemes considered previously 
indisputable (between them, all aspects of the Boltzmann description even for the week- 
disorder limit in the treatement of the low temperature resistivity, 1 ) have been forced to be 
re-analyzed. Based on these facts, suggestions to follow new roads have been made, under- 
lining that the disordered materials cannot be understood by evading the real issue, and 
forcing the disorder into a mold of procedures standard for ordered systems 1 . Furthermore, 
it has been stressed that the non-perturbative view on disorder could lead to significant ad- 
vancement in the understanding and description of such syatems 10 . In the same time, several 
recent developments in the field require the non-approximated solution of the wave-equation 
for disordered and interacting systems as a key feature for a much deeper understanding of 
the emerging processes, and their interpretation, especially where it is expected, or exper- 
imentally is seen, that the electrons are maintaining their long-range phase coherence and 
retain their wave nature, as in the case of solid grains, short wires, dots, mesoscopics, prox- 
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imity to critical points, presence of long-range order, or of some kind of order in general 2 ' 13 , 
presence of quantum interferences 14 , etc. 

On this background, the first steps towards exact results for disordered systems have 
been made. On this line especially non-periodic models of different type were analyzed. 
In these models, the non-periodicity is considered as introducing the effect of the disorder 
in the Schrodinger equation, and is taken into account in different ways, for example as 
non- analytic behavior in the potential 15 , incommensurate potential 16 , quasiperiodicity 17 , 
topological disorder connected to tesselation 18 , local bond-orientational order 19 , etc, these 
possibilities presenting also interdependences between them. This way of describing the 
disorder cannot be considered as representing the level of simple toy models only, since 
besides the fact that real physical systems holding such properties are known 2,15,16,19 , there 
are concrete cases where it is also known that a such type of representation (for example 
through quasiperiodicity) give analogous behavior for the system as random or disordered 
potentials 15,16 . 

In D — 1 the majority of studies leading to exact results were given for Fibanocci 
type of lattices 15,20,21 , and the interested reader will find more extended information on 
the ID subject in review articles, like Ref. 2 . For D > 1, which is for interest in this 
paper, the first exact results have been obtained for quasicrystalline systems, where first, 
theoremes dealing with structure have been formulated, such as those involving inflation 
rules, or Conway's theorem 22,23 . Later on, in few cases, even exact eigenfunctions have 
been deduced for Penrose lattice 22 in 2D 24 ~ 27 . This type of lattice being a prototype of 
quasicrystalline systems 28 , clearly exceedes the level of a curiosity of pure mathematical 
character since it is related to nearest-neighbor bond orientational order which is observed 
for example in simulation of supercooled liquids and metallic glasses 19 , attracting clear 
interest 29 . In the 2D Penrose lattice, in a simplified view, flat and thin rhombuses cover the 
plane completely, forcing the resulting pattern to be non-periodic and introducing disorder 
in the system. For these systems, in 2D exact eigenstates were obtained by Kohmoto and 
Sutherland 24 for a strictly localized state including in the Hamiltonian on-site disordered 
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potential as well depending on the number of bonds entering in a given site, Sutherland 25 
obtaining a self-similar state taking into account as well on-site potential which may have 
eight different values depending on the nature of the site, Arai et al. 26 obtaining new strictly 
localized states in comparison with those described in 24 , and Repetovitcz et al. 27 obtaining 
eigenstates even in the presence of plaquette-diagonal hoppings. The knowledge of these 
results has clarified puzzles related to the influence of the disorder in several aspects, in an 
extent which only exact results can provide. We note on this line clarifications of disputes 
originating from the interpretation of numerical results 24 , evidences for the self-similarity of 
some ground-states in disordered systems 25 , evidences for singular features of ground-states 
in certain non-periodical systems 26 , occurence of degeneracy proportional to the system 
size in eigenstates 26 , scaling properties of the exact ground-states 25 , relative stability of the 
confined states on boundary conditions 26 , possible existence of allowed and forbidden sites 
in the eigenstates 26 , etc. 

We underline that the above mentioned exact eigenstates are valid only for independent 
electrons, e.g. they were deduced from models build up on a tight-binding Hamiltonian in 
r space describing a single electron moving on an aperiodic graph 17 . Because of this reason, 
and in light of the facts previously presented, it would be extremely stimulating for the field 
to see in what extent the deduced exact properties at independent electron level, remain or 
not valid in exact terms for a really multielectronic and interacting system as well in the 
presence of disorder. In our knowledge, exact results of this type, up to this moment, are 
not known. 

In this paper, we report for the first time exact ground-states depending on the electron 
concentration for multielectronic and interacting 2D systems in the presence of disorder. 
The ground-states are paramagnetic, lose their eigenstate nature in the independent elec- 
tron approximation, present properties known in the Penrose-lattice (for example strong 
degeneracy proportional to the system size), describe a localization-delocalization transition 
driven by the electron concentration, and in the localized case present clear evidence for 
long-range phase coherence. 
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The results are reported for two-band type models. The presence of two bands is not di- 
minishing the applicability of the results since, from one side, real materials are of multiband 
type, and the theoretical description is given usually by projecting the multiband structure 
in a few band picture 30 , which is stoped only for its relative simplicity at the one-band 
extreme level, when this is possible. From the other side, the experimental one, several ma- 
terials treated traditionally in a two-band picture have been experimentally found containing 
disorder and presenting extremely interesting properties (as non-Fermi liquid behavior for 
example 31 ) whose emergence is considered connected to the presence of the disoreder (see 
for example 32,33 and cited literature therein). Concerning the presence of real random sys- 
tems holding two type of electrons, we mention the intense activity related to rare-earth and 
actinide compounds which behave as random Kondo insulators 34 ' 35 ' 32 holding two type of 
electrons (d and /), whose properties are described in a fixed lattice topology, but randomly 
distributed Hamiltonian parameters 33 . 

The procedure we use originates from developments leading to the first exact ground- 
states for the periodic Anderson type models obtained at finite value of the interaction in 
ID 36-38 , 2D 39 ~ 41 , and 3D 42 , which have been made here applicable in the disordered case as 
well. Our model is build up on a 2D graph in r space, whose all vertices are of the same 
rank (four edges are collected by every vertex), so the topology is fixed. Four neighboring 
nearest-neighbor vertices form elementary plaquettes, and hopping (including the non-local 
hybridization as well) is possible along the edges and diagonals of elementary plaquettes. 
On each vertex local on-site potentials are acting, and on each vertex local on-site Hubbard 
type repulsion is present as interaction. 

The remaining part of the paper is build up as follows. Section II. presents the Hamil- 
tonian, Section III. describes an exact transformation of the Hamiltonian which allows the 
deduction of the presented results, Sect. IV analyzes the disorder present in the system and 
provides concrete examples for the emergence of the model conditions necessary for the so- 
lutions to occur, Sect.V. presents the exact ground-states, Sect. VI describes ground-state 
expectation values, and finally, Sect. VII. concluding the paper closes the presentation. 
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II. THE HAMILTON! AN OF THE MODEL. 



The fixed topology of the described system allows us to treat the problem in a 2D tight- 
binding Hamiltonian defined in r space on a 2D Bravais lattice with disordered Hamiltonian 
parameters. For this system we consider an unit cell / described by the primitive vectors 
(x, y), and we take into account two type of electrons denoted by the particle index p as 
p = d, f. In these conditions our starting Hamiltonian has the form H = H + Hm, where 



^o= E EE 

p=d,fp'=d,f v 



r^O 

h**= e E (i) 

P=d,f i 



In the one-particle part H , the ,, length" of the hopping denoted by r with possible non- 
zero values x, y, y + x, y — x, is allowed to extend only to distances contained in /, i.e. 
nearest-neighbors (x, y) and next-nearest-neighbors (y + x, y — x) (see Fig.l). Denoting 
by iVA the number of lattice sites in the system, te random nature of H is given by a) 
the 2A^a independent, non-correlated, random (repulsive) on-site Hubbard interactions U[, 
p — d, f contained in Hi nt , and b) 2N^ new independent, non-correlated and random H 
parameters chosen (as will be clarified below) from the (site, direction and spin dependent) 
t\f +r r cr amplitudes. We underline that the t p,p ' coefficients contain hybridization (p ^ p'), 
and on-site potential (r = 0) terms as well. 

We demonstrate below, that for certain local conditions imposed for t p ' p ' parameters 
which maintain the number of AN\ independent non-correlated random variables in the 
system, the exact multielectronic ground-state wave function of H in the interacting case 
can be explicitly given in an electron concentration dependent manner. 

We mention that the spin-dependent nature of the Ho parameters is not essential for our 
deduction. The H parameters can be in principle spin dependent as well, and we underline 
this aspect, in order to extend the applicability of the results also to Hamiltonians with 
non-diagonal hoppings 43 ' 46 ' 47 too. Furthermore, concerning the type of the model we use, we 
mention that for U- 1 = 0, the H from Eq.(l) represents a disordered periodic Anderson model 
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(or Anderson lattice), while for ^ 0, Eq.(l) describe a disordered two-band Hubbard 
model. Our results are applicable in both cases. For physical realization of such type of 
systems see for example Ref. 34 . 

We further consider that the mobility of the two type of electrons present in the system 
(d and /) is different, and the ratio in mobility is the same on all lattice sites. As a 
consequence, from the point of view of hopping amplitudes, starting from amplitudes written 
for d electrons, we have 

t?f' = «/i>./ + V,/t?> d (2) 

fc i,i+r,r,<r w fc i,i+r,r,(r J ) 

where w is a (site independent) measure of the mobility ratios between / and d electrons. 
We mention that hopping amplitudes between different orbitals often satisfies such type of 
relations in real systems 47 . 

Concerning again the t v,v ' terms, being interested in the behavior of particles given by the 
disordered hoppings and interactions, we only consider situations for which the localization 
of particles in local trapping centers is avoided, i.e. we have 

*S5v, > o • (3) 

In the following Section we are presenting a transformation of H in a form that allows to 
obtain exact ground-states in its spectrum. 

III. THE TRANSFORMATION OF THE HAMILTON! AN. 

Let us introduce a numbering of the lattice sites by the integer number / in the studied 
2D lattice containing N A = L x L lattice sites, strarting from the down-left corner in the 
lowest row (I = 1), going from left to right up to the end of the firts row (I = L), then 
going upward and continuing with the second row again from left to right, and so on. In 
this manner, for example around an arbitrary lattice site i, numbered by / = i, we find 
the site numbering notations presented in Fig. 2. a. The introduced notation allow us to 
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turn from a vectorial site notation to a scalar one, which simplify as well the notation of 
the Hamiltonian parameters. For example the t?f +rra for r = x, (r = y) at site I = i 
becomes tff +lx a ., {t\f +Ly (T ) respectively (see Fig. 2. a.). Similarly, the next-nearest neighbor 
components (x + y, y - x) become ^f +1+L>x+Vi<T , t^ i+Lna . 

Let us further introduce a plaquette operator A- i a defined for every arbitrary cell 1^ taken 
at site i (see Fig.2.b). The cell I[ is denoted by its down-left corner i. The sites inside Ji, 
are numbered in a cell independent manner by the index n — 1, 2, 3, 4 starting from the site 
i and counting anti-clockwise inside the unit cell I\ (see Fig.2.b). In these conditions we 
obtain for A- ha the expression 

4 

A\,a = ^2{a n ,dd i+rnt0 . + a n j/i +rn)<T ) , (4) 

n=l 

where a ntP are numerical coefficients, the same in all unit cells (detailed description of this 
procedure can be found in Ref. 40 . For unit cell independent notation of the coefficients a n>p 
see Ref. 42 ). Let further connect to every unit cell I[ =i , two random variables e^, and e^. 

Our results are based on the observation that if we define the plaquette operator param- 
eters a n p via the non-linear system of equations 

td,d _ * i * 

i,i+l,x,a — a \ : d a ^,d e i,c ~r a4 t d a 3,d e i-L,a , 

td,d _ * i * 

i,i+L,y,a ~ a l,d a i,d e i^ + a 2,d a 3,d e i-l,a , 

,d,d _ * 

l i,i+l+L,x+y,<j ~ a l,d a 3,d e i,v , 

,d,d _ * 

^i+l,i+L,y-x,a ~ a 2,d a 4,d e i,CT ) 

±d,d _ | |2 , I 1 2 , I 1 2 , I 1 2 /r\ 

l i,i,0,<r — \ a l,d\ e i,v + \ a 2,d\ + |03,d| Q-1-L,ct + \0>A,d\ ei-L,cr , [0 ) 

and a„j = wa n ^ holds, where the parameter w (see Eq.(2)) is real but arbitrary, then, 
taking into account periodic boundary conditions, the one-particle part Hq of the starting 
Hamiltonian from Eq.(l) becomes 

H = J2^tA^- (6) 

i,cr 
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Compairing the last equality of Eq.(5) to Eq.(3), we obtain the condition e ij(T > 0, although 
e^ a are random variables. As a consequence, H in Eq.(l) becomes positive semidefinite 

# = E + E E W, A , (7) 

i,<r P=d,f i 

and this property preserves the potential possibility to obtain the explicit form of the ground- 
state in the interacting case. 

IV. THE DISORDER IN THE SYSTEM. 
A. The presence of randomness in the model 

Before going further, we should analyze the kind of randomness we have in the system. 
We start with the observation that H in Eq.(7) which will be analyzed further on, is clearly 
disordered since contains 4Aa independent, non-correlated (non- negative) arbitrary random 
variables ej iCr and Uf. However, the randomness must be understood not only at the level 
of the transformed Hamiltonian Eq.(7), but also at the level of the starting H presented in 
Eq.(l). Since the disorder in H int is the same in Eq.(l) and Eq.(7), this question relates 
only the randomness in H . In order to understand the source of the disorder in H from 
Eq.(l), we have two different alternatives. 

One possibility for this, is to observe the linear relationship between the on-site energy 
levels tff Q(J and e^ a in the last row of Eq.(5). As a consequence, we can consider that the 
initial disordered parameters of the starting H in Eq.(l) are the rj it<T = t^f (T variables whose 
number is and the e ijtJ new disordered parameters from Eq.(7) are obtained from these 
by a linear transformation 

_ | 1 2 , I \2 , I |2 , I |2 (q\ 

Vi,cr — \0'l,d\ e i,cr + |«2,rfj Q-l,cr + \ a 3,d\ e i~l-L,a + \0-A,d\ ^i-L,a , [°) 

which contains also 2N\ equations. Since the r] itU disordered parameters are independent, 
in this view one can consider that H, besides the randomness in Hint, possesses as well 
,, diagonal-disorder" at the level of H in Eq.(l). 
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Alternatively, one can consider in Eq.(5) the unit-cell diagonal hopping amplitudes 
(tf,i+L,y-x,(Ti tif+i+L,x+y,<j) directly proportional to e iiCr as the source of the disorder in the 
one-particle part of the Hamiltonian, H Q in Eq.(l). In this case, H is considered to contain 
besides the randomness in Hint, also ,, non-diagonal" disorder at the level H Q . 

In both cases, the remaining equalities in Eq.(5) must be considered as local constraints 
necessary for the solutions to occur. Since the number of H parameters in Eq.(l) is much 
higher than the number 2N\ of random one-particle variables, these constraints do not 
alter the random nature of the disordered variables ({Uf , rji ta }, or {Uf,e it(7 }). Rather, they 
lead to I.) interdependences between H parameters not containing the disordered variables, 
and II.) connect other H parameters to rj ij0 - or e ij(T disordered variables. These constraints 
emerge in the process of the transformation of Eq.(l) into Eq.(7), and we underline that our 
solutions are valid only in the case when this transformation can be done (i.e. Eq.(5) holds). 
Both cases mentioned above as non-diagonal and diagonal disorder in H will be analyzed 
in details below. 



B. Connections to the solutions obtained for Penrose tiling. 

Considering the disorder in H as non-diagonal, the here presented solutions can be 
viewed as arising from extension of the conditions used in the exact study of the Penrose 
tiling. In order to understand this statement, let us introduce the constants K\ = a\ d a^ d , 
K 2 = o,2,d a 4:,dy an d observe, that since e iiCr are random, the diagonal (next-nearest neigh- 
bor) hopping matrix elements in every unit cell U, namely tf* = tfi +1+Lx+ya . and 
tf~ = ti+ li+Ly _ x 0, excepting their ratio (Ki/K 2 ), remain random as well 

tf,a — Ki£i,<r ? tile — K 2 e it<T . (9) 

Considering for example the hopping amplitudes without directional dependence, i.e. tf 2 a = 
tile, and taking for simplicity K 1 = K 2 = 1, we obtain 

<2,. = ^ CT , (10) 
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which (excepting the fixed sign of e^ a > 0) means completely random and independent unit- 
cell diagonal hoppings for all spins in all unit cells (see Fig. 3.). As a consequence, based on 
Eq.(9) or its particular form from Eq.(10), we see that the randomness given by {Uf,ei tU } 
in Eq.(7), can be considered to originate from the randomness given by {Uf,tf^} at the 
level of the starting H presented in Eq.(l). In this case, once the hopping amplitudes along 
the diagonals of every unit cell have been randomly chosen, the remaining t p,p ' parameters 
can be determined based on them. The study of Eq.(5) shows that fixing the tf^ values, we 
have the liberty to choose independently two more constants relating the one-particle part 
of the Hamiltonian H , namely K 3 , and K4, (signK > 0, K = K2K3K4), based on which 

,d,d Ty- ,d,d 1 1 ,d,d ,d,d _ jy- ,d,d . ^ ■ d,d 

^i,i+l,x,a ~ t^&i-L+l^^-x^cr ' jy- ^i,i+l+L,x+y,a^i,i+L,y,a ~ 4^i,i+L-l,y-x,a ' ^i,i+l+L,x+y,a i 

,d,d _^2_.d,d ^_.d,d ^U,d,d K"i_.d,d 

^i,i,0,a j£ ^i,i+L+Y,x+y,a ' j£ ^i—L—l,i,x+y,a ' j£ ^i,i+L—l,y—x,cj ' j£ ^i— L+l,i,y— x,a i V / 

and the numerical coefficients present in Eq.(5) in function of K m , (m = 1,2,3,4), ar- 
bitrary parameters become a\ d a 3 ^ — Ki, a 2 d a 4,d a l,d a ^d — 
K1/K4, a*2,d^,d = K 2 K 4 , al d a 34 = K 2 K 3 , |a M | 2 = Kf/K, \a 2 , d \ 2 = K 2 K 4 /K 3 , \a 34 \ 2 = 
K, \o-A,d\ 2 = K 2 K S /K 4 . In order to have real value for all t p ' p ' parameters, all K m must be 
real. To understand in details Eq.(ll), let us introduce short notations as well for nearest- 
neighbor and local amplitudes in the form t*£ = tjiw *?,'<? = tfi+L^,*, *f° = *mJd,<t» which 
represent the tf% amplitudes for d electrons with spin a in unit cell Jj for v = |r|. Using 
these notations, Eq.(ll) becomes 

t% = K 3 tt~ L , a + K^tf; , tgf = K 4 tt\, + K^tf; , 

where R\ = K\jK and R 2 = K^jK 3 . As shown in Fig.4., the amplitudes presented 
in Eq.(12) are determined by the t d,± unit cell diagonal amplitudes that surround tf%. For 
example, as seen from Fig. 4. a, the tf* nearest neighbor hopping amplitude (full line arrow) 
is given by the tf„ and tf^L a random unit cell diagonal amplitudes (dotted line arrows) 
which start from the same site i and intercalate tf%- Similar situation is present for tf% 

11 



in y direction (see Fig.4.b), while tf® as seen in Fig.4.c. is determined by the four t d,± 
,,plaquette-diagonal" amplitudes that start from the same site i. 

Concerning Eq.(12), we mention that in the study of disordered systems, constraints 
(correlations) between bond and site properties are often considered. The constraints a 
priori introduced can be in some cases even of long range type, as taken for example in the 
case of isotropically correlated random potentials 48 , correlated networks 49 , etc., and even 
calculation techniques have been developed in order to deal with ,,constrainted" disorder, 
for example in the form of correlated random numbers 50 , or random matrices with symmetry 
properties or holding constraints 49 . In our case, local constraints exist which connect the 
plaquette diagonal bond hoppings (considered as the true independent random variables of 
H Q ), to edge (nearest-neighbor) bond hoppings and local one-particle potentials. Since the 
plaquette diagonal bond can be unambiguously connected to the plaquette, in the described 
case, random plaquette properties (i.e. random bonds connected to plaquettes), through 
local constraints presented in Eq.(12), fix nearest-neighbor or local amplitudes. 

Conrete physical situations where in disordered systems the random plaquette properties 
determine nearest-neighbor or local amplitudes are also known in the literature. For exam- 
ple, in the case of topologically disordered system of Caer type 44 using random mosaics, 
very similar to Voronoi tessellation generated from disordered arrangement of particles 51 , 
random flips of plaquette-diagonal bonds performed with a given probability, determine the 
local nearest-neighbor hoppings, and introduce in this way the disorder in the system 18 . 
Concerning disordered on-site one-particle terms generated by random bonds connected to 
plaquette properties, we mention for example the Penrose lattice 22,23 case, where the on-site 
one-particle potentials have been introduced by the local coordination number 24 . In the 
mentioned case, practically the random on-site one-particle potential at site i is determined 
by the number of bonds entering in the site i. Our on-site potential 1%£ Qa given in Eq.(12) 
and presented in Fig.4.c. has clear similarities with this choice, since reduces to a such 
type of behavior in the case in which all tf^ unit cell diagonal hoppings are equal, and, for 
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example, for all m, K m = 1. The diference between Eq.(12) and Ref. 24 from the point of 
view of the random on-site potential is that in our case, the on-site potential is determined 
by the value of the random bonds entering in the site, while in 24 , by the number of the bonds 
entering in the site. So contrary to Ref. 24-26 , where the study has been concentrated on the 
effects of the lattice topology alone, in this paper we analyse the problem in a fixed topology, 
concentrating on random t p ' p ' values. We further mention, that in the Penrose lattice case, 
when also the plaquette diagonal hopping amplitudes are taken into account at the level 
of exact independent electron eigenstates 27 , solutions are found only when constraints are 
present between hopping amplitudes. 

Let us consider a concrete physical example in support of Eq.(12) which demontrates as 
well that solutions deduced in the context of Penrose lattice 27 use quite similar conditions. 
For this, let us take a simple spin-independent case t P ' p ra = t P ' p r _ a = t P ' p r , and consider 
a situation for which, randomly positioned A or B atoms in the middle of the elementary 
plaquettes providing the random unit cell diagonal hoppings tf^ = tf ± introduce the ran- 
domness in H Q . In this situation, tf^ is either t d,± (A) or t d,± {B), depending on the type 
of atom situated in the middle of an unit cell. For this example Eq.(12) express the fact 
that the hopping amplitude along a bond (nearest-neighbor hoppings tf x and tf y ) depends 
on the randomly situated atoms A and B placed in the neighborhood of the bond, and 
that the on-site energy of a given site (tf °) depends on the randomly positioned A and B 
atoms in the neighborhood of the site, which are physically quite acceptable conditions. The 
linearity of these interdependences can be physically motivated by the small influence of the 
atoms A or B not situated directly on the bond or on the site. For the presented exam- 
ple, t d,± (A) (t d,± (B)) in the Penrose lattice case would correspond to the notations di,d 2 , 
(g^,^) used in Ref.( 27 ). Furthermore, in Ref.( 27 ) tf x = tf v = 1 is considered, and our t d '° 
is denoted by Cj. The conditions in which solutions are obtained for the Penrose lattice case 
(see Eqs.(3.9),(3.11),(3.13) in Ref.( 27 )) are in fact linear relations of the type of our Eq.(12). 
The main difference between our model and that of Ref.( 27 ) at the level of H is that in 
our case, the plaquettes described by t d,± (A) and t d,± (B) can emerge completely random, 
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while in Ref.( 27 ), the „plaquettes" („rhombi") described by (di,d 2 ) and (d 3 ,d 4 ) emerge only 
with the randomness allowed by the Penrose tiling. Because of this reason, our solutions 
extend the exact solution possibilities known in Penrose lattice case to non-quasicrystalline 
disordered systems even in the presence of electron-electron interaction. 

C. The disorder seen in the one-particle part of H as diagonal-disorder. 

Considering the source of the disorder in H diagonal, the random parameters of the 
model become Uf and r] i (T . In this case Eq.(5) requires two supplementary conditions to be 
satisfied as follows. I.) The next nearest-neighbor hoppings surounding a nearest-neighbor 
hopping (see Fig.4.a,b, and first two equalities from Eq.(12)) must be related at the level of 
hopping amplitudes. These conditions are not specific for the presented disordered model, 
but are rather connected to the method itself. Indeed, such conditions we find solving non- 
disordered cases as well (see Ref.( 39 ~ 42 )), and the obtained hopping amplitude ratios are 
delimitating parameter space regions where the obtained solutions are valid. II.) The next- 
nearest-neighbor (unit cell diagonal) hoppings starting from a given site i are all together 
linearly related to the on-site energy level (considered disordered here) at the site i (see the 
last equality from Eq.(12), or alternatively Eq.(8), and Fig.4.c). This local constraint, in 
this form, is specific for the random case studied here. 

Even if the conditions I. and II. presented above seems to be quite specific at first view, 
we show below that they ar compatible with the presence of the diagonal disorder on physical 
grounds. For this to be visible, we analyze a simple pedagogical spin-independent hopping 
case which is x-y symmetric as well, so t£* = tf£ = fj£ =1 = tf\ tf^ = tf^ a = tj£ =2 = tf 2 , 
and K 3 = K 4 = K , R 2 = 1 (see Eq.(12)). Let further consider for the study that the random 
on-site potential rji is created by the randomly positioned A T atom at site i of the lattice 
with lattice spacing a, where the index r fixes the type of the atom. In this manner, if the 
atom A T will be placed at site i, it creates the on-site energy level rji = r\ T . After this step, 
we must model the expression of the distance dependent hopping amplitude U(r) for the 
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electron which starts the hopping from i. Taking into account a simple exponential distance 
decrease, we may simply take tj(r) = CiBi(Av)e~ ar , where the constant a describes the 
distance decrease (r ^ is considered). The amplitude CiBi(Av) is build up from the 
component Cj = C T depending on the energy level at site i (the atom A T present at site i), 
and the average effect of all surounding atoms felt at site % denoted by B^Av). Since only 
tf 1 and tf 2 hoppings are considered, we have for the tf 1 , (tf 2 ) case the r = a, (r = a\/2) 
argument value in U(r). 

After these considerations, the conditions I. and II. mentioned above look as follows. 
Condition I. links together 3 hopping amplitudes for hoppings which start from the same 
site % (see Fig.4.a,b, and the equalities relating tf 1 from Eq.(12)), providing the condition 
e aa(V2-i) = j{ + (1/Ko). As can be seen, condition I. determines in fact the strength of the 
hopping (parameter a) through the constant Kq, and is not specific for the random 
mentioned above. Rather, it fixes the t d,1 /t d ' 2 ratio introducing limits for the valability of 
the solutions in the T = phase diagram of the starting Hamiltonian. As a consequence, 
we can further analyze condition II. considering a known parameter. 

Condition II. (Fig.4.c, and the equality relating rji = tf® from Eq.(12)) links together 
four next-nearest-neighbor hoppings which again start from the same site. As a consequence, 
taking into account that the atom A T is placed on the lattice site i, we find i] T = C T Bi(Av) (2+ 
i?i + (l/ J R 1 ))e- v ^ aa . The remaining Bi(Av) coefficients must be deduced at each site from 
the condition i* ■ = tjj. As can be seen, condition II., through the parameter R\ relates the 
disordered on-site energy level values to the hopping amplitude components Cj. 

As presented above, the conditions necessary to be fulfilled for the solutions to emerge are 
present in disordered systems, being compatible to a truly acceptable physical background. 
Taking into account more complicated parametrizations for tj(r), the equation of r\i^ = tf® 
in Eq.(12) (since has in its right side the same r value), reduces to an equation for the 
amplitudes of the r function in U(r), while the remaining equalities in Eq.(12) determine 
the ti(a)/ti(ay/2) ratios. 
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V. THE GROUND-STATE WAVE-FUNCTION. 



Starting from the positive semidefinite structure of H in Eq.(7), the ground-state wave 
function \fy g ) is obtained for H\^ g ) = 0. Now, let us concentrate first to the H component of 
H presented in Eq.(6). Taking into account Eq.( 4), and as shown under Eq.(5), a n j = wa nt d, 
where w is arbitrary but real, we realize that 

4 

Ai,a = £ a n,d( d i+r„,cr + wf i+Tnt<T ) , (13) 
n=l 

so in the right side of Hq in Eq.(6) only operators of the form Oj^ a = (d^ a + wfj j0 ) are 
present. If now we define 

Ol = d],« - , (14) 

which satisfies Oj >a Oj l(jl = —0j, ,O^ a ^ then taking = lljt^^a + v jdj,-o\\ty, where Hj 
is taken over different (although arbitrary) lattice sites, Vi are arbitrary coefficients, and |0) 
is the bare vacuum with no fermions present, we obtain 

£ ei , ff itAj*> = • (15) 

i,cr 

Since introduce fermions (d or /) with arbitrary spin, strictly on different sites, double 
occupancy is avoided, and 

E £^AI*> = o, (i6) 

P=d,f i 

holds as well. Since the minimum possible eigenvalue of H in Eq.(7) is zero, the ground- 
state for arbitrary N < N^, where iV represents the number of electrons within the system, 
becomes 

w=n[C+ w AUii°>- ( 17 ) 

In Eq.(17), the product must be taken over j sites which can be arbitrary chosen, and 
different j values must be related to strictly different lattice sites. The ground-state wave 
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function of H given in Eq.(7) for A < A A (i.e. at and below quarter filling) can be always 

written in the form of Eq. (17). As a consequence, for A = A A the ground-state becomes 

Na ~-t ~-t 

\* g (N = A A )> = H[O ja + VjO^WO) . (18) 
i=i 

For A < N A , since the j sites in Eq.(17) can be arbitrary chosen, the complete ground-state 
become 

\* 9 (n < n a )) = j2{a RN n tC + v A,-«\}\°) . ( 19 ) 

where the sum J2r n is made over all different domains containing A < N A lattice 
sites from the system, and olr n are numerical coefficients. Furthermore, it is important to 
underline that Eqs. (18,19) represent the ground-state only in the interacting case (at least 
one of on-site two-particle interactions Uf must be non-zero at all sites i, since otherwise, 
because of the presence of the double occupancy, \^ g ) in Eq.(17) is no more an eigenstate 
of the Hamiltonian) . 

In my knowledge, Eqs. (18, 19) contain the first exact multielectronic ground-state wave- 
functions obtained in 2D for a disordered system in the interacting case. As explained above, 
Eqs. (18,19) are no more eigenstates for the independent electron approximation. Since the 
ground-state in the interacting case, even for infinitesimal interaction, changes qualitatively 
in comparison to the non-interacting case, Eqs. (18,19) cannot be connected in a perturbative 
way to the ground-state of the disordered but non-interacting system. 

The GSs presented above are strongly degenerated. Their degeneracy at quarter filling 
is given by the N arbitrary Vi values and the arbitrary (but non-zero) w, is proportional 
to the size of the system. The existence of such type of states for 2D Penrose type lattices 
has been first conjectured by Semba and Ninomiya 53 and Kohmoto and Sutherland 24 , and 
further analyzed in 26 ' 54 . From the here reported results it can be seen that this property is 
present also for other systems as well in the multielectronic and interacting case too, at least 
for N = N A - We stress however, that in the case N < N A , the degree of the degeneracy 
strongly increases given as well by the geometrical degeneracy present in Eq.(19). The order 
of magnitude of the degeneration becomes in this case N R = N A \/[N\(N A — A)!]. 
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VI. GROUND-STATE EXPECTATION VALUES 



A. The localized case. 

Despite the possibility to chose the Hamiltonian parameters in a spin-dependent way, the 
obtained GS is globally paramagnetic. At 1/4 filling (N = N A ), the GS contains rigorously 
one electron on each site, so the hopping is completely forbidden in GS, and as a consequence, 
the system is localized, holding long-range density-density correlations. 

Indeed, calculating the ground-state expectation values through Eq.(18), in this case we 
find for arbitrary i ^ j and all a, a' 

(4,4.') = 0, {f}Jis) = 0, (d{J^} = 0, {f}Jjs)=0, (20) 

where (...) = (^ g (N A )\...\^ g (N A ))/{^ g (N A )\^ g (N A )), \^ g (N A )) is presented in Eq.(18), and 
(^ g (N A )\^ g (N A )) = n£\[(l + M~ 2 )(l + \v%\ 2 )]- Th e reason for Eq.(20) is simple: \V g (N A )) 
contains exactly one electron on each lattice site, so \^i{p,p')) = p ij(T Pj a j\ i ifg(N A )), where 
p,p' = d, f and i ^ j, contains a double occupancy, and as a consequence \^i(p,p')) and 
\^J g (N A )) are orthogonal. We underline that since Eq.(20) holds for arbitrary it remain the 
same after the average over the disorder (i>j variables). Denoting the translational invariant 
averages by ((...)) = / P({vi})(Y\ i dv i ) (...), where P({t>;}) describes the distribution of the 
disordered variables (being here arbitrary) and / P({vi})(J\idvi) = 1 holds by definition, 
Eq.(20) automatically implies as well {(p\ a p'j, a :)) = for all p,p' = d,f, all a, a' and all 

Furthermore, introducing for % ^ j the notation D(i,j) = [(l + \w\~ 2 )(l + \vi\ 2 )(l + \vj\ 2 )] 2 : 
we find 

I r^j . J 2 I _ j 2 

AiA°) = H^D(i,j), (ra£_ ff n£_ ff ) = ' J D(i,j), 

W-A,-«) = ^^(^ K-A,«) = \^ D ^j)- (2i) 
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Starting from Eq.(21), for hi = ^(raf^ +nf jCr ), based on Eq.(21) one obtains 

(hihj) = 1, ((n^-)) = 1, (22) 

where the second equality holds as explained below Eq.(20). Introducing now Sf = 
(1/2) [(ra^ + nf T ) — (hf^ + nf^)], based again on Eq.(21), for % ^ j we have 

(5?5/> = ^ (23) 

where p n = (1 — |f n | 2 )/(l + \v n \ 2 ) takes arbitrary values in the domain (— 1,+1), so 
{{S*Sj)) = arises. As can be observed, \^ g (N\)) indeed describes a paramagnetic, com- 
pletely localized ground-state containing long-range density-density correlations. At quarter 
filling, since \^ g (N = N\)) coherently controls the particle number occupancy at all lattice 
sites forbidding the hopping (and non-local hybridization) in the same time, the GS clearly 
presents phase coherence over the whole lattice. 

B. The delocalized case. 

Under quarter filling, empty sites emerge in the GS, Eq.(20) deduced through Eq.(19) 
does not hold, hopping is no more forbidden, and as a consequence, a derealization oc- 
curs, the system becoming itinerant (remaining further paramagnetic). Indeed, in this 
case, at A < N\, based on Eq.(19), the GS wave function can be written as \^ g (N)) = 
Er n <* Rn \*(Rn)), where \*(R N )) = |* Rjv (M)> = ({««}) IW^L + Vjd[_ a )\0) 

build up an ortho-normalized wave function set containing Nr components, we have 

jei?jv(l + \ v j\ 2 )i an d ®r n are coefficients independent on the 
disordered {vi} set. The operators of the type p\ a p^ a , where p,p' = d, f, now have non-zero 
matrix elements between ground-state components \^(Rn)), \^(R' n )) describing R N , R' N do- 
mains of the form R^ = -Djv-i + i, R' N = -Dat-i + j, where D^-i represents an arbitrary 
region of the lattice containing A — 1 lattice sites, and i, j are representing two different but 
arbitrary (not necessarily nearest-neighbor) sites of the lattice. We have for example 
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V(i + kl 2 )(i + N 2 )' 

(*(^)|/U>I*»)> = 7= (1 , + ,'T)?" 1 | ^ • 

V(i + I^r)(i + I^l 2 ) 

Since the disordered variables emerge in Eq.(24) through \vi\, \vj\ non-negative numbers, the 
average over the disorder maintains the non-zero values in Eq.(24). As a consequence, the 
hopping being no more forbidden, the system becomes indeed itinerant. Since as seen from 
Eq.(24) all d or / electrons can hop everywhere in the ground-state, the wave function in 
Eq.(19) is clearly an extended state. The conducting nature of the extended states can be 
demonstrated (see for example [ 55 ]) through the variation of the chemical potential as the 
number of electrons vary. For this reason, let us observe that the ground-state wave function 
from Eq.(19) acting on the Hamiltonian from Eq.(7), by the construction of the wave function 
as explained in Sect.V., has the property H\V g (N < N A )) = E g (N)\^ g (N < N A )) = 0, 
where E g (N) is the ground-state energy for iV particles in the system. Since \^ g (N < N\)) 
is a wave function with non-zero norm, this relation means E g (N) = 0. As a consequence, 
for TV < N\ — 1, we have for fj,+ = E g (N + 1) - E g (N) and pT = E g (N) - E g (N - 1), the 
expression 

^+ - = 0. (25) 

Therfore, the state we analyze is conducting (see also [ 56 ]). 

Furthermore, the SfSj operator will have non-zero matrix elements only along the di- 
agonal in i?7v variables, and in conditions mentioned for pi after Eq.(23) we further have 
«5f5/»=0. 

As can be seen, Eqs. (18,19) describe a localization-delocalization transition driven by the 
electron concentration p n , which emerge at p c n = 1/4, the delocalized phase being present 
in the region p n < p c n . The occurrence of this transition is intimately connected to the 
multielectronic nature of the description which is made in the presence of the inter-particle 
interaction and absence of trapping centers. Indeed, the problem considered at the level of 
independent electron approximation (e.g. absence of inter-electronic interaction) in presence 
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of trapping centers leads to a one-particle problem in the presence of an attractive potential, 
which ends up usually at small energies with localization. Here all these are avoided. 

Concerning the possibility of the emergence of Griffiths phases in influencing the de- 
scribed transition, we mention that the Griffiths singularities arise due to the presence of 
statistically rare clusters that are anomalously strongly coupled, and hence they are unique 
features of the disordered system (see for example [ 57 ]). The effect becomes weaker with 
increasing dimension, increasing interaction, increasing the number of the components N 
of the dynamical variables. In our case N = 3 (for example for spin), the dimension of 
the (quantum mechanical) description is D = 2, and the results are not valid at zero inter- 
electronic interaction. All these conditions make unlikely the major influence of Griffiths 
phases, especially when the Pi parameters are all maintained perfectly random at all sites 
as mentioned below Eq.(23), prohibiting in this way the local formation of anomalously 
strongly coupled clusters. 

VII. SUMMARY AND CONCLUSIONS. 

We deduced exact multielectronic concentration dependent ground-states for disordered 
and interacting two-dimensional quantum mechanical systems at and below quarter filling. 
The ground-states describe a localization-delocalization transition driven by concentration 
and provide paramagnetic behavior. The ground-state nature is lost in the absence of 
the interaction e.g. independent electron approximation. The deduced results are non- 
perturbative and cannot be perturbatively reached from the non-interacting, altough disor- 
dered case. The studied system is of two band type, and the disorder is present independently 
in both Hi n t and H parts of the Hamiltonian, the trapping centers being excluded. The pre- 
sented procedure extends the exact solution possibilities known in 2D Penrose lattice case to 
non-quasicrystalline disordered systems as well, even in the presence of the inter-electronic 
interactions. 
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FIGURES 

i+y i+x+y 




i+x 



FIG. 1. Unit cell I at site i. The vectors x, y denote the primitive vectors of the unit cell, and 
arrows indicate the possible r values allowed for the hoppings. 
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i+L-1 i+L i+L+1 



a 



i-L-l i-L i-L+1 



b) 




i+x+y 



i l+X 

FIG. 2. (a) The numbering of sites in a L x L two dimensional lattice around the lattice site 
i, and (b) The unit cell I{ placed at an arbitrary lattice site i, together with the i independent 
notation (n = 1, 2, 3, 4) of sites inside I\. (x, y) are the primitive vectors of the unit cell. 
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b) 
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/+ 







t d,0 

s 

Ii-L 



i-L-1 i-L 



FIG. 4. The amplitudes at site i for a) v = x, b) v = y, and c) v = respectively. is 
denoted by full arrow in a),b), and by a circle in c). Dotted arrows with ± label indicate the tj'^ 
random amplitudes that enter in the expression of tf^ presented in the plots a),b),c) from Eq.(12). 
In all plots the notation of the unit cell Ij at site j, containing t^ v a is also presented. For example, 
in a) Ii defined at site i contains if'j", and Ii-L defined at the site i — L contains tf~ T . 
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